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We describe here the properties expected of a higher (with emphasis on the order fourth) order 
phase transition. The order is identified in the sense first noted by Ehrenfest, namely in terms 
of the temperature dependence of the ordered state free energy near the phase boundary. We 
have derived an equation for the phase boundary in terms of the discontinuities in thermodynamic 
observables, developed a Ginzburg-Landau free energy and studied the thermodynamic and magnetic 
properties. We also discuss the current status of experiments on BaosKoABiOs and other BiOz 
based superconductors, the expectations for parameters and examine alternative explanations of the 
experimental results. 



I. INTRODUCTION 

There has been a proposalii^ recently that the transi- 
tion to superconductivity in the cubic, non-copper based, 
moderate Tc (32K) material Bao.eKoABiOsiBKBO) is 
of order four in the sense proposed by Ehrenfest"^. The 
free energy and its first three derivatives are continuous 
across the phase transition. The fourth derivatives with 
respect to temperature (or, alternatively, with respect to 
some mechanical variable such as pressure or magnetic 
field) are discontinuous. This is somewhat unexpected. 
After nearly a century of an extensive study of phase 
transitions, where the transitions were either discontinu- 
ous (in the Landau sense of a discontinuous change in the 
order parameter and/or entropy, the latter in the sense 
Ehrenfest envisioned) and therefore first order, or con- 
tinuous where the transition was believed to be II order, 
a consensus seems to have emerged that all departures 
from a strictly II order behavior can always be explained 
in terms of thermodynamic fluctuations. 

To recall, the Ehrenfest definition of the order of a 
phase transition can be summarized in the thermody- 
namic free energy F{B,T). Here B is the magnetic field 
and T is the temperature. Thus near the phase boundary 
(for a superconductor) the free energy satisfies: 



F(B,T)« [l-T/T,] 



[1 - B/B,2] 



(1) 



Here p is an integer, the order of the transition and 
fi and C are small corrections. The scaling field Bc2 ~ 
a^^, (a = ao{l — T/Tc)) defines the exponent X2- Thus all 
derivatives of this free energy, of order less than p, will 
be continuous, the p*'' derivative will be discontinuous if 
/i and C are zero. Otherwise they will be divergent with 
exponents /i for a thermal derivative and X2C for a field 
derivative. For a second order phase transition, p = 2. 
The lowest thermal derivative, which may be singular, is 
the specific heat (or the magnetic susceptibility). Apart 
from conventional superconductors and fermionic super- 
fiuids where mean field theory is valid and fi and C are 
both nearly zero, most transitions are accompanied by a 
small but finite /i. 

That fi is small, is on the one hand a validation of 



the Ehrenfest classification. On the other hand, it fa- 
cilitates formalism, such as the renormalization group, 
to calculate the numerical value of fi and Ci depending 
on symmetry properties alone. Disorder brings about 
an uncertainty to this clean picture. Imagine for exam- 
ple a system with a distribution of T^s caused by spa- 
tial inhomogenieties. The goal here is to establish that 
a broad transition may be more than just a dirt effect. 
The transition may also be higher order in that the ex- 
ponent p above is larger than 2. It might be tempting 
to use the known formalism for a second order transition 
to calculate larger (in magnitude, since they are nega- 
tive) values of a fi. However, since the mean field theory 
for a higher order transition is different, the calculations 
should appropriately be carried out within the framework 
of a higher order transition. 

In a second order phase transition, the slope of the 
phase boundary in the B-T plane is related to the dis- 
continuity in specific heat (AC) and in the magnetic sus- 
ceptibility (Ax), 



AC - T,Ax 



dB. 



c2 



dT 



(2) 



If AC is zero then either the phase boundary has a zero 
slope or Ax is also zero and the phase transition is of 
higher order. There are materials where AC = at the 
superconducting transition. 

The class of cubic perovskite materials, the bis- 
muthates, seem to be the promising candidates for the 
effects discussed here. Following the discovery.^, of super- 
conductivity at IIK in BaPbxBii-x03(BPB0), there 
has been a continuing interest in this class of materi- 
als. While BKBO keeps the record^ for the highest Tc 
(32K) and the Rubidium-doped samples Bai-xRbxBiO^ 
(BRBO, rc=25K)'^'® have T^'s which are shghtly smaller, 
most of the materials in this series have a transition tem- 
perature nearer to lOK. They are all missing the specific 
heat anomalies. 

The new materials in this class are K and Rb doped 
SrBiOs with = 12 - 13i^2, and Ki^xBh+xOs with 
a = lOK^". The straggler in this group is BP SO 
{Ba{Pb, Sb)03)ii. with the highest of 3.5K, but other- 
wise with similar properties as the corresponding ones in 
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the other members of the group. These materials await 
further work such as a careful measurement of specific 
heat and magnetic properties as a function of magnetic 
field and pressure. 

The natural question here is what happens if the two 
discontinuities are simply very small. This is similar to 
what happens in a second order phase transition which 
turns weakly first order. The entropy should have a kink, 
instead it undergoes a small jump. Here, too, the system 
approaches with all the accoutrement (for example the 
temperature and field dependence of the free energy as 
noted above in Eq.(^) of a higher order transition, but 
then turns weakly second order. The suggestion below 
is that the order of the transition should be identified 
from the temperature dependence of the free energy as 
one approaches the transition. 

This paper is a follow-up to an earlier brief letter^ on 
the report of and an interpretation of the anomalous su- 
perconducting transition in BKBO. It is also a theoreti- 
cal complement to a detailed report of the experiments^. 
The conclusions of ref.(0) are based on three independent 
results in the magnetic field and temperature dependent 
magnetization of BKBO. They are (recall that a here 
refers to the reduced temperature (1 — T/Tc)): 

1. The thermodynamic critical field has the tem- 
perature dependence (all fields are in Tesla) Bo{T) = 
0.51ai-8i. 

2. The lower critical field has the temperature de- 
pendence with a cross-over. Near Tc, it is Bci{T) = 
0.096a^ °^. At low T, the temperature dependence is lin- 
ear in a. 

3. The upper critical field Bc2{T) = ig.Ta^ ^i. The 
Ginzburg-Landau parameter k = \jBc2/Bci = 14.4a°'^^. 

These results, as we discuss below, challenge the view 
that the phase transition to superconducting state is of 
second order. They instead lead to an internally consis- 
tent description of a fourth order phase transition. Re- 
sult (1), in the context of Eq. shows that p = 4 and 
fjL = 0.4. Result (2) is a consistent with the field theory 
developed for a fourth order phase transition. In general, 
we can show that for a p*'' order transition, London pen- 



etration depth (and therefore Hd) satisfies A « a 



(p-i). 

We can also show that within mean field theory, the su- 
perconducting coherence length satisfies f ~^ w a. Thus 
the parameter k = A/^ is constant only for a II order 
phase transition. In a transition of order higher than 2, 
K diverges at Tc- Physically, one might view this as a 
separation between the magnitudes of flux penetration 
and order parameter variation. Near Tc, flux expulsion 
gets weaker for a higher order transtion. Althought as 
we see below in Sec. (IV), the estimates derived for low 
temperatures correspond to k « 14, a large value. 

The paper is organized as follows. In sec. II, we dis- 
cuss the thermodynamic foundations of a fourth order 
phase transition. Sec. Ill contains discussion of a model 
field theory for a higher order phase transition. Sec. IV 
follows up with a discussion of the gradient terms in the 
field theory under the heading: magnetic field effects. 



Sec. V is a discussion of the reported results and the 
existing estimates of the specific heat discontuity. While 
they seem to agree, there are reasons to believe that the 
agreement might be conincidental. Sec. VI concludes 
with a summary of principal results and a discussion of 
assumptions. 



II. THERMODYNAMICS OF A HIGHER 
ORDER TRANSITION 

In the Ehrenfest definition of the order of a phase tran- 
sition, p^^ order phase transition is characterized by a dis- 
continuity in the p*^ derivative of the free energy. Con- 
sider the free energy as a function of the temperature T 
and a mechanical variable, say the magnetic field B. The 
discontinuities also determine the phase boundary in the 
B — T plane. When p=l, the discontinuous derivatives 
are the entropy S and the magnetization M . We have the 
Clausius-Clapeyron equation: dB/dT = - AS /AM. If 
the entropy and the magnetization are continuous across 
the phase boundary, the transition may be of higher or- 
der. For p=2, [dB/dTY = AC/T^Ax. 

It is possible to generalize these results. For a third or- 
der phase transition, the corresponding relation becomes 



dB 



dT 



AdC/dT 
T.AdxIdB 



(3) 



To derive this relation, note that in the vicinity of the 
phase boundary separating phases 1 and 2; I'Hopital's 
rule will necessarily change Eq. (01 in to 



dB 



dT 



AdC/dT 
T,Adx/dT 



(4) 



Since the susceptibility x is now continuous across the 
transition, it follows that: 



dB 
dT 



^dx/dT 
Adx/dB 



(5) 



Multiplying Eqs.Q) and ©, we recover Eq.Q. 

Similarly, the expression for a fourth order phase tran- 
sition becomes: 



dB 
dT 



Ad'^C/dT^ 
TcAd^x/dB"^ 



(6) 



A generalization to an arbitrary order phase transition 
is nearly self-evident here. The proof is by induction, just 
as we derive Eq.® from © or Eq.Q from Eq.©. 



dB 
'dT 



AdP-'^C/dTP-' 
TcAdP-^x/dBP- 



(7) 



The expression is not unique, although certain of its 
aspects are notable such that (a) the slope is raised to 
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the power of the order of the transition and (b) that 
the numerator consists of the discontinuity in the ther- 
mal derivative of the free energy while the denomina- 
tor contains the derivative with respect to the mechan- 
ical variable, here the magnetic field. Of course, there 
are other expressions involving mixed thermal and me- 
chanical derivatives which have lower powers of the slope 
dB/dT. There are also relations between mixed deriva- 
tives corresponding to the well known Maxwell relations. 
The expressions above seem to be the most compact and 
symmetric. 

It is possible to write the corresponding expressions for 
the case when the mechanical variable is pressure. For 
a second order phase transition, the density remains the 
same across the transition but the compressibility, the 
second derivative of the free energy as a function of pres- 
sure is discontinuous. Since the longitudinal sound veloc- 
ity directly involves the compressibility, if the transition 
temperature depends on the pressure, the sound veloc- 
ity changes at the transition. In a fourth order phase 
transition, the compressibility is continuous across the 
transition as well as its first derivative, but its second 
derivative, with respect to both the temperature or the 
pressure is discontinuous and should be observable in the 
temperature/pressure dependence of the sound velocity.. 

In the above, the implicit assumption has been that fi 
and ^ are zero. If they are nonzero, then the discontinu- 
ities in free energy derivatives turn into weak logarithmic 
divergences. One can then take the exponents for these 
divergences and derive scaling laws between them. As- 
suming that the magnetic field and temperature depen- 
dences in a p*'' order phase transition are described by 
Eq. J^l, Kumar and Saxena^^ have derived the scaling 
laws. Defining the exponents as: 



{p-l)li + p[3 + K = p{p-l) (8) 
K = P[{p - 1)5 - 1] (9) 
l^L = p-P[5 + 1] (10) 

{5+l)n+[{p-l)5-l]{^J,-p) = Q (11) 

This brings us to the end of the section on thermodynam- 
ics. The next section develops possible alternatives for a 
field theory, underlying the thermodynamic properties. 



III. A MODEL FIELD THEORY 

We need to construct a free energy F{il)}, a function of 
the complex order parameter ip = Ae''^ such that when 
F is minimized with respect to the order parameter ip^ 
its minimum is the thermodynamic free energy. F is also 
a measure of the fluctuations, i.e. those configurations of 
ipir) which are near the equilibriunm configuration i}}o{r). 
In the spirit of Ginzburg-Landau formalism, F could be 
written as an expansion in the powers (even powers, since 
F must be gauge invariant; a change in </) to </) -f c, where 
c is a constant, should leave F invariant) of ij). However 
such an expansion necessarily leads us to a second order 
phase transition. 

To consider a different temperature dependnece of 
the free energy, suppose we write the leading term in 
F « a?/'^* (recall here that a — ao{l — T jT^ and t is de- 
fined here as an exponent), then to be consistent with Eq. 

we need to make sure that = a*^P^^^ This can 
be accomplished in seemingly two different ways. First 
consider the choice, ■0^ ~ Qi then t = (p-1). The free en- 
ergy for a fourth order phase transition, in this approach 
can be written as 

F{i) = Ae*"^} = -aA^ -f 6A« -I- Fq (12) 

where Fq represents the terms of various powers in 
jV'i/'p. The second term with the power 8 ensures the 
first of the conditions, that IV'P = a. 

The second possibility is provided by the choice t = 1, 
leading to a free energy (to avoid confusion, here we use 
n for the order parameter), 

F{n:\^~ar? ^hri^l^ ^Foin} (13) 

Eqs. (|12|) and H13|l are related via the transformation 
n = tp^ . An ambiguity here is the choice of the cube, for 
example, the transformation could well be n = iV'PV'j 
then the phases for -0 and n are identical. We have 
no way of deciding on this choice. The reader is left 
to speculatei^ about the possibility of a fractional flux 
quantum. 

So far we have not considered the gradient terms. In 
Eq. H13|l . it is clear that it must be c|Vnp. Transformed 
back into Eq. H12|l . the gradient term becomes |V'0^p. 
Even in Ea. H12() it can be argued that the gradient term 
must be as described here. If it were not and say the 
gradient term were simply |V?Ap, the temperature de- 
pendence of this term would be a^. In the presence of 
a magnetic field the transition would switch to a second 
order. Assuming that « a, the gradient terms consis- 
tent with a fourth order transition are, 

c,\y^p'\' + c2\v^^^\^ + c3\y'M^ (14) 

The additional terms are indicative of a nonlinear 
Meissner effect and will be omitted from any further dis- 
cussion of stability effects. 
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We now have a free energy which can also be written 



as 



F{^{r),T}= dVij'^[-a^p^ + b^* + c\Vi:\^] (15) 



This looks like the usual free energy, multiplied by 
an external weight factor proportional to "0^. In the 
BCSi^ii^ derivation of the Ginzburg-Landau free energy, 
the prefactor is the density of states, which determines 
the energy scale of the phenomena. Here the extra factor 
may be seen as arising from a density of states which is 
proportional to ip"^. Thus the effective density of states 
in the normal state vanishes. It is non-zero only in the 
condensed state where it has a given temperature depen- 
dence related to the order of the phase transition. In 
effect we have a transition from an insulator to a super- 
conductor. 

Szabo et al— have measured the temperature depen- 
dence of the energy gap by point-contact tunneling with 
a silver electrode. They find a BCS like temperature 
dependence. This would be consistent with the descrip- 
tion above if A were to be identified as the energy gap 
(or more precisely, it rules out the possibility of n being 
identified as an energy gap). 

Let us also note that the free energy for a third order 
phase transition is given by (corresponding to Eqs. H12|) 
and ifT^ . 



-I- bn^ + c| VnP 



(16) 



In the above equations, the relation is n = t/;^. 

Finally let us note that the interaction term with the 
power 8/3 in Eq. I|13|l may well arise from an averag- 
ing carried out over some other field. In the spirit of 
a Hubbard-Stratonovich transformation, there may be 
auxiliary fields which interact with the superconducting 
order parameter and are ultimately responsible for the 
unusual order of the transition. Looking at Ea. (|12|) . one 
might get the impression that it is necessary to find some 
physical reason for the absence of terms quadratic and 
quartic in ip. The answer is to some extent, provided by 
Eqs. p3ll5|l . The primary message here is to point out 
that there are many different ways to see the mathemat- 
ical structure necessary for a higher order phase transi- 
tion. 

In a real physical system, there are often nearby in- 
stabilities. For example in BaBiOs, there are several 
charge-density-wave (CDW)iii^ transitions (see for ex- 
ample Ref.(@)). The corresponding transition tempera- 
tures decrease with K concentration, x and the highest 
of them vanishes near x=0.4 . This value of x coincides 
with the highest superconducting Tc- A microscopic the- 
ory may well involve a subtle interplay between the CDW 
and superconducting degrees of freedom. 

To conclude this section, note the results arising from 
free energy minimization of Ea. (|13|l . 



3a 

4&' 



for T<Tc 



(17) 



The thermodynamic free energy (the free energy at the 
minimum) Fo{T) is given by; 



Fo{T) 



27 a-* 
250 63 



Fm{T) 



(18) 



AC - C{T) - Cn{T) = -T-^iFo - F^) 



AC 



81 a 



64 Trb^ 



^{l-T/T,f for T<T, 



(19) 



The specific heat approaches Tc with a power law. 
Here for a fourth order transition, the power is quadratic. 
There is a discontinuity in d'^C/dT'^ at Tc; 



A 



d^C 81 at 



dT^ 32 r3&3 



IV. MAGNETIC FIELD EFFECTS 

This section addresses the consequences of the gradi- 
ent terms in the free energy in Eas. (|12ll3|l . In a charged 
superconductor, the effect of magnetic field appears via 
the gauge transformation, V — > V -I- 2TriA/(po, where A 
is the vector potential and ipo = h/2e is the supercon- 
ducting flux quantum. The free energy Eq. lfT5|l is also 
supplemented by the field energy density term, _B^/2/i. 
Focussing on issues specific to this section, note that the 
London equation is the Euler-Lagrange equation for the 
vector potential A{r). Following Eq. 112() . 







dF 
dA 



2^0 



-V'A- 



,27r, 



'AM = 



(20) 



The flux is expelled whenever A ^ 0. There is a cor- 
responding penetration length, A such that (recall that a 
is the reduced temperature, a » (1 — T /Tc)) 



X-^T) 



2/ioc(27r/(^o)^ A^ = a3 



(21) 



This is a mean field result for a fourth order phase 
transition. For a general order p, of the phase transition. 



the temperature dependence of A 



j(p-i). This is 



quite distinct from the fluctuation contributions to the 
exponent xi (defined as A^^ « a^^). For a second order 
transition, the classical value of a;i = 1 but the fluctu- 
ations tend to make it smaller. For example, for a 3d 
XY model, the superfluid density exponent xi « 2/3, 
measuredi^ in several high Tc superconductors. For a 
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p order transition, fluctuations will probably reduce it 
to < p — 1 but not by much and the connection to the 
classical value and the corresponding order, as well as the 
contrast with 3d XY model should remain clear. 
The lower critical field-'^^ in a superconductor Bci{T) 



lip 1) This should be a clear observable. In- 



deed the expression for Bd (T) contains corrections loga- 
rithmic in K = A/^ where ^ is the temperature dependent 
coherence length. In a second order phase transition, 
both A and f have the same temperature dependence 
making k a constant. In a higher order transition, n 
is temperature dependent and thus the temperature de- 
pendent exponent for Bci{T) has logarithmic corrections; 
even at the mean-field level, it may be slightly smaller 
than {p— 1). 

The divergence of k indicates a novel physical effect. 
As the condensation occurs, the flux expulsion takes 
place more gradually than the development of order pa- 
rameter stiffness. There is no reason that they should 
happen together (as they do when the transition is of 
second order). Indeed it could be argued that the curi- 
ous behavior is the constancy of k. The range of values 
of K in BKBO is such that the system is always a su- 
perconductor of type II, even at low temperatures. This 
may also explain why the measured values of k vary so 
much between different experiments. They have been 
measured at different temperatures. 

Let us note that for a third order transition, the tem- 



perature dependence of A 



a^. If the transition in 



Bi-2212 is of third order^Sli then the penetration depth 
will correspond to a different power law. The exponent 
may be smaller than 2 due to the fluctuations, but it 
would be considerably larger than one, the number char- 
acteristic of a second order phase transition. For a gas 
of bosons, Schafroth^^ has derived a temperature depen- 
dence for A^^ « a^lna. There are conflicting reports 
of two experiments, one noted-- the temperature depen- 
dence as quadratic in the reduced temperature, the other 
noted^^ its absence. 

Based on the experimental results^ for BKBO, which 
include the lower critical field, Bci(r) = 0.0955a^ "^ 
Tesla, the upper critical field Bc2 = 19. 7a^'^^ Tesla and 
the thermodynamic critical field Bo — 0.51a^'^^ Tesla, we 
can derive the characteristic limiting length scales at T 
= 0. 



A(0) 

m- 



t^M - 830A 



(22) 



Finally, in a second order phase transition supercon- 
ductor, there is an identity B^ = BciBc2 ■ This identity 
is not satisfied here. The geometric mean of the two criti- 
cal fields (in Tesla) is 1.37a^'^^ while the thermodynamic 
field, as noted above is O.Sla^'^^. There is no clear rea- 
son to draw any conclusions from this observation. It 
may be that the identity is simply not expected to be 
valid in general. 



Finally, let us note that the Euler-Lagrange equation 
for the order parameter is simply derived. 



= -6atp + Sbi// - 4 



2VV (23) 



We will defer a detailed study of the consequences of this 
equation to until later. 



V. DISCUSSION OF EXPERIMENTS 

In retrospect, it is interesting to note that soon after 
the original discovery- of superconductivity in BPBO by 
Sleight et al, Methfessel et alM' published a paper with 
the title "Why is there no bulk specific heat anomaly at 
the superconducting transition temperature of BPBO" . 
The authors first wondered about the nature of super- 
conductivity in BPBO, whether it was some sort of im- 
purity phase since magnetic measurements of supercon- 
ducting transition have often turned up to be misleading. 
But they also noted that "if it is a bulk effect, then we 
may have to invoke a new, as yet unidentified, mecha- 
nism for its cause" . Soon after the discovery of super- 
conductivity in BKBO by Cava et al§., the early specific 
heat measurement3SS*S& were already reporting the miss- 
ing specific heat anomaly. There was one reporliSi of a 
small specific heat discontinuity but the two other mea- 
surements were unable to verify it. 

It is possible to estimate the size of the expected dis- 
continuity in specific heat. When this question is put 
before the cognoscenti, the answer one is more likely to 
get is that the specific heat discontinuity in BKBO is 
small, of the order of a few mJ/moleK. This estimate, as 
we discuss below, is based on questionable assumptions. 
Most of the reports of measurements of this magnitude 
are less than convincing. 

We show below that the estimate is based on arguable 
evidence. For instance, several authors-* have used for- 
mulae based on thermodynamic identities for a second 
order phase transition, in one case, a version of the fa- 
miliar Ehrenfest expression, Eq. (O where the identity 
Ax = [Sttk^]"^ converts this relation into an expression 
for AC in terms of experimentally known quantities k 
and the slope of the Bc2 curve near Tc- Kwok et al2£ used 
this expression to estimate AC/Tc w 3.75 mJ/mole/K^. 
In the measurements of Hall et aU*2i, k apparently di- 
verges near Tc as k w (1 — t)^^. If we use this tempera- 
ture dependence, we get AC = 0. 

Similarly, based on a more fundamental expression. 



AC = r, 



92 



Bl 



(24) 



where Bo is the thermodynamic critical field. Batlogg et 
alii estimated AC/T^ « 2.1 ± 0.4 mJ/mole/K^. Once 
again, if we use the experimentalljiii^ measured temper- 
ature dependence of Bo « (1 — T/TcY, we get AC = 0. 
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One begins to get the impression that the small value of 
AC may well be an imperfect estimate of zero. 

On the other side of this story, there is an estimate 
which yields an expected AC at least three orders of 
magnitude larger. According to Hundley et al2S*2&, the 
normal state specific heat in BKBO can be fit to an elec- 
tronic contribution of C/T = 7 = 0.15 J/mole K^. The 
data, when plotted as C/Tys-T"^, show two straight lines 
with a crossover around 16K. On the low temperature 
side, the T=0 intercept is zero. However on the high 
temperature side (> 16K^ transcending Tc), the inter- 
cept has the value quoted above. If we take the view 
that the low temperature side is the propertji^S of the 
superconducting state, then the effective 7 derived from 
the high temperature side must be the large value. Since 
AC w jTc, we expect a discontinuity at the supercon- 
ducting transition of the order of 4-6 J/moleK. 

In a recent letter, Woodfield et al2i have questioned 
all of the conclusions of ref. (0) on the grounds that they 
(Woodfield et al) have measured a 7 = Q. 9m J/moleK^ 
and AC/Tc = 2.1mJ/moleK^, a result in excellent 
agreement with the estimates of Batlogg et al. But, as 
the discussion above suggests, there may be a different 
point of view for these small estimates for AC. In any 
case, a non- vanishing discontinuity in a higher derivative 
of AC still leads to a structure below Tc- For a fourth 
order phase transition, the discontinuity in second ther- 
mal derivative of the specific heat means that the specific 
heat rises below Tc, although quadratically in the reduced 
temperature a. It eventually vanishes at T = 0. That 
means that there must be a specific heat maximum below 
the transition temperature. It is just that the maximum 
(or some measure related to it, such as the temperature 
half way to it) should not be interpreted as an estimate of 
Tc- Nor is the height of the maximum to be interpreted 
as an estimate of the discontinuity. Indeed, whereas in a 
second order phase transition, the size of the discontinu- 
ity is a condensate property, the corresponding informa- 
tion in a fourth order transition lies in the coefficient of 
the term in the specific heat. 

There is a report from the Grenoble group'^^ of mea- 
surements of specific heat and the magnetization in 
BKBO. This report emphasizes the presence of a specific 
heat maximum below Tc, the latter derived from tun- 
neling measurements associated with a vanishing energy 
gap. Near Tc, the specific heat appears to be quadratic 
in a. Moreover, they have obtained a measure of 
from the magnetization which seems to fit a^'^. If the 
transition is second order, the mean field exponent is 
one. Fluctuations, within a 3d XY model yieldi^ smaller 
exponent (of order 2/3). This result is also in sharp con- 
trast to the results of HalU*2i, where the exponent is closer 
to 3. This discrepancy, amongst experiments, of a fac- 
tor of 2 in the exponent, is unexpected and needs to be 
addressed. 

There are reports of energy gap measurements from 
tunneling experimentsiSi^ which satisfy the BCS char- 
acteristics. Thus the ratio 2A/kTc is in the range of 



4 — 4.3, only slightly larger than the BCS value (3.4) and 
possibly indicative of strong coupling effects. There is 
no reason to expect higher order transition effects in the 
temperature dependence of the order parameter (which 
is presumably also the energy gap). For example the 
temperature dependence of the energy gap measured by 
Szabo et aiiSi very well fits BCS expectations. Ea. (|15l 
describes an order parameter which is conventional in its 
temperature depndence, but the free energy temperature 
dependence is not. 



VI. SUMMARY AND CONCLUSIONS 

To recall, perhaps the more convenient definition of 
the order of a transition is the one based on Eq. 
namely it is the integer p, in the exponent [p — /i in the 
temperature dependence of the free energy. This leads to 
the following conclusions. 

(1) The superconducting phase transition in BKBO 
has the characteristics of a IV order phase transition in 
the sense described by Ehrenfest. This indentification is 
based on the temperature dependence of the themody- 
namc critical field, the temperature dependence of the 
lower critical field (or the London penetration length) 
and is consistent with the missing discontinuity in mag- 
netic susceptibility and the specific heat. 

(2) The equation for the phase boundary in ap*''-order 
phase transtion is described by Eq. Q . The correspond- 
ing relation for a 4*''- order transition, Eq. 10) relates 
the lowest order singular derivatives of the free energy. 

(3) A model field theory, consisting of polynomials in 
the order parameter magnitude, is able to describe the 
anomalous temperature dependences in item (1) above. 
In particular, the free energy leads to the anomalous tem- 
perature dependence of the thermodynamic as well as 
the lower critical field. This free energy may arise in the 
vicinity of competing phase transitions 

The possibility of a higher order phase transition be- 
comes viable if the specific heat discontiuity is absent. In 
BKBO^^, as well as in BRBO^ (Rubidium instead of 
Potassium) and in BPBO (Pb replacing Bi)^**, the spe- 
cific heat discontinuity is anomalous, either missing or 
believed to be extremely small. Several discussions in 
the literature advocate the belief that, in BKBO, AC 
is indeed expected to be quite small. We have analyzed 
these estimates and find that the small estimates of AC 
may well be imperfect estimates of zero. When combined 
with the other indications (anomalous temperature de- 
pendence of the thermodynamic critical field and of the 
lower critical field as well as the missing discontinuity in 
the magnetic susceptibility) the evidence for an unusual 
transition is notable. This evidence is strictly thermody- 
namic. The line of reasoning is model- independent. 

It is possible to understand this transition in terms of 
a model free energy whose expansion in powers of the 
order parameter begins with the 6*'* order term. If we, 
as is commonly done in textbooks, assume that all phase 
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transition must be described by a Ginzburg-Landau free 
energy with an expansion in all powers of the order pa- 
rameter (including quadratic and quartic), then this as- 
sertion is equivalent to saying that higher order phase 
transitions do not exist. That the thermodynamic evi- 
dence described here is flawed in some unknown but fun- 
damental way and must be rejected. If on the other hand, 
the experimental evidence for anomalous temperature de- 
pendence stands the test of further experiments, then we 
may have to abandon the "all powers" requirement for 
the free energy. One possibility might be gleaned from 
Eq. H15|l . Here the free energy density is an "all powers" 
expansion, although the total free energy is an integral 
over all space with a "weight" function. 

There remains the task of identifying the microscopic 
origins of the phenomena discussed here. At this stage 



the experimental evidence does not provide clear insight 
into an answer to this question. Certainly, future experi- 
ments will. Hopefully the above discussion will stimulate 
more experiments in future. 
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